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Abstract 

A new oriented method using a diazonium salt reaction was developed for linking /^-adrenoceptor (/J 2 -AR) on the surface of 
macroporous silica gel. Stationary phase containing the immobilised receptor was used to investigate the interaction 
between /? 2 -AR and ephedrine plus pseudoephedrine by zonal elution. The isotherms of the two drugs best fit the Langmuir 
model. Only one type of binding site was found for ephedrine and pseudoephedrine targeting /i 2 -AR. At 37 °C, the 
association constants during the binding were (5.94±0.05)x10 3 /M for ephedrine and (3.80±0.02) x10 3 /M for 
pseudoephedrine, with the binding sites of (8.92±0.06) x10~ 4 M. Thermodynamic studies showed that the binding of 
the two compounds to /i 2 -AR was a spontaneous reaction with exothermal processes. The AG 9 , AH 6 and AS 9 for the 
interaction between ephedrine and jS 2 -AR were -(22.33±0.04) kJ/mol, -(6.51 ±0.69) kJ/mol and 50.94±0.31 J/mol-K, 
respectively. For the binding of pseudoephedrine to the receptor, these values were — (21.17±0.02) kJ/mol, — (7.48±0.56) 
kJ/mol and 44.13±0.01 J/mol-K. Electrostatic interaction proved to be the driving force during the binding of the two drugs 
to /? 2 -AR. The proposed immobilised method will have great potential for attaching protein to solid substrates and realizing 
the interactions between proteins and drugs. 
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Introduction 

The interactions between a protein and a drug are of growing 
interest due to the role in elucidating protein functions, explaining 
drug action mechanisms and discovering novel drug candidates 
[1,2]. G-protein coupled receptors (GPCRs) are the largest and 
most diverse superfamily of integral membrane receptors due to 
their crucial role in curing varied diseases [3-5] and being the 
most important class of drug targets [6] . 

Currently, a wide range of methods have been developed for 
investigating the drug-protein interaction, including ultrafiltration 
[7], mass spectrometry [8], X-ray crystallography [9], nuclear 
magnetic resonance [10], surface plasmon resonance [11], affinity 
capillary electrophoresis [12] and computer-aided methods [13]. 
Although contributing gready to the exploration of drug-protein 
interactions, these approaches need to be improved to overcome 
the limitations of long analysis time, requirements for samples with 
high concentrations or purities, potential interfering substance 
from the sample and relatively specialized equipment [14]. High 
performance affinity chromatography (HPAC) is another widely 



used method for the study of drug-protein interaction [15]. The 
method often involves immobilization of a protein on a solid 
matrix to construct affinity stationary phase. 

Synthesising an assay for attaching a protein to a solid matrix is 
a key factor for accomplishing an HPAC method. Physical 
absorption and random immobilised methods are commonly 
reported ways to synthesis an affinity stationary phase [16,17]. The 
types of stationary phases constructed by these methods have been 
confirmed to generate the inevitable issue of losing bioactive 
binding sites of the protein during the immobilisation procedure. 
Oriented immobilisation has proved feasible to address this issue 
in previous publications [18,19] and should be further used to 
synthesise HPAC stationary phases attributed to the particularly 
important role in life sciences [20,21]. 

Beta 2 -adrenoceptor (/? 2 -AR), one of the members of the 
GPCRs, is a crucial target of many drugs for fighting obesity 
and ailments of the heart and respiratory system [22]. Ephedrine 
and pseudoephedrine (Figure 1) are isomers widely used for 
treating diseases of the respiratory system. Both ephedrine and 
pseudoephedrine exert their therapeutic action through the same 
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signal pathway involving /? 2 -AR [23]. The two drugs have varied 
potency in vivo although they have similar structures and the same 
acting mechanism. One reason for this may be the different 
binding affinities of the two ligands to /? 2 -AR. This work is firsdy 
designed to establish an oriented immobilised method for fixing 
his-tagged fusion /J 2 -AR on a solid matrix to synthesise an HP AC 
stationary phase. Further work is performed to investigate the 
binding behaviours of ephedrine and pseudoephedrine to /3 2 -AR 
by HPAC methods. 

Theory 

Zonal elution 

The binding of ephedrine and pseudoephedrine to the oriented 
immobilised j8 2 -AR was determined by zonal elution. The theory 
assumes that a known concentration of a competing agent (I) 
continuously passes through a column containing an immobilised 
ligand (L) while the injecting volume of an analyte (A) is small 
enough to be neglected compared with the concentration of I in 
the mobile phase. Assuming I and A have direct competition at a 
single binding site on L and the interactions have association/ 
dissociation kinetics, Eq. (1) can be used to calculate the binding of 
A and I to L [24,25]. 



1 
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(1) 



where X is the capacity factor for A due to the existence of a 
second type of binding site on the column, and K is the overall 
capacity factors of A ascribed to all the types of sites. In this case, 
the equation assumes the capacity factors of A produced by the 
individual sites are additive, and an independent estimate of the 
capacity factors due to the second type of site can be obtained, 
allowing for the capacity factor attributed to the first type of site 1 
(Site 1) to be measured by k{ = /c'-X. 

When the column is equilibrated with an eluent at concentra- 
tion [A] and [Tj equal to zero, and the retention of compound A is 
injected as a perturbation, Eq. (2) is used to determine the binding 
of A to L: 
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Assuming a single type of binding site, a linear expression can be 
obtained by plotting 1 / sqrt (k') as a function of [A] to determine 
m L and K A . For two types of binding sites, a second term is added 
to the k' expressions. 
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Figure 1. Chemical structures of ephedrine and pseudoephed- 
rine. (a), /fi,25-ephedrine; (b), 7S,2S-pseudoephedrine. 
doi:1 0.1 371 /journal.pone.0094955.g001 



Thermodynamic study 

The thermodynamic interaction between /? 2 -AR and ephedrine 
plus pseudoephedrine was described by the van't Hoff equation, 
listed as Eq. (3) and (4). 



\nK A - 
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In these equations, R is the ideal gas constant (8.31 J/mol/K), 
and T is the temperature in Kelvin. According to the semi- 
empirical rules summarised by Ross et al. [26], when AG e <0, 
AH e >0 and AS 8 >0, hydrophobic interaction is the main 
interaction force during the binding process of a drug to protein. 
If AG 9 <0, AH e <0 and AS 6 >0, electrostatic interactions mainly 
contribute to the interaction. If AG e <0, AH 6 <0 and AS 9 <0, Van 
der Waals forces or hydrogen bonds are the major forces of the 
binding. 

Materials and Methods 

Chemicals and instruments 

Ephedrine (li?,2^-ephedrine, 171241) and pseudoephedrine 
(7>S,.2.S'-pseudoephedrme, 171237) were obtained from the National 
Institutes for Food and Drug Control (Beijing, China). Fatty acid 
and globulin free bovine serum albumin (BSA, A0281) was 
obtained as lyophilized powder from Sigma-Aldrich (St. Louis, 
MO, USA). The Quaternary Sepharose Fast Flow pre-packed 
column (16x25 mm), Ni 2+ Sepharose 6 Fast Flow (10x58 mm) 
and Sephadex G25 pre-packed column (10x100 mm) were 
purchased from GE Healthcare Life Sciences (Uppsala, Sweden). 
Macroporous silica gel (SPS 300-7, pore size 300 A, particle size 
7.0 mm) was from Fuji Silysia Chemical Company Limited 
(Tokyo, Japan). 

The HPLC system employed in this work consisted of a binary 
pump, a column oven from Agilent Technologies (Santa Clara, 
USA) and a diode array detector (Waldbronn, Germany). The 
data were collected and processed by a Chemstation 5.2 software 
installation. 

Expression and purification of /? 2 -AR 

The full gene of Sus scrofa /? 2 -AR (Refseq:MN_ 001128436) 
was cloned in our previous work [27] and the protein was 
expressed in E.coli BL21 (DE3) with an N-terminal His-tag. The 
transformed E. coli strain was grown overnight at 37 °C in 50 mL 
Luria-Bertani (LB) medium containing 100 Hg/L penicillin. When 
the value of OD 60{) reached 0.4—0.6, isopropyl-/3-D-thiogalacto- 
pyranoside was added to a final concentration of 2.0 mM, and the 
culture was incubated at 30°C for 6.0 h. The cells containing His- 
tagged /? 2 -AR were separated from the culture media by 
centrifugation at 6000 xg for 30 min at 4.0 °C, and the 
supernatant was discarded. Falcon tubes containing the cell pellets 
were placed into a dry ice/ ethanol bath for 3-5 min and were 
transferred to cold water for 60 s to melt the ice along the walls of 
the tubes. The pellet was broken free by a sterile plastic spatula 
and suspended in 200 mL neutral lysis buffer (50 mM Tris-HCl, 
pH 7.4, containing 0.3 M NaCl and 0.5 mM EDTA) in a 500-mL 
flask. Lysozyme (0.5 mg/mL) was added to the buffer and 
incubated for 2.0 h on ice with agitation to digest the cell wall. 
The E. Coli cells were lysed using 20 mL of a 10 x Mg/Mn salt 



PLOS ONE | www.plosone.org 



2 



April 2014 | Volume 9 | Issue 4 | e94955 



Binding of Ephedrine and Pseudoephedrine to fo-AR 



solution (100 mil MgCl 2 and 10 mM MnCl 2 ). DNA was removed 
by DNase (10 (ig/mL) with incubation for 30 min at RT. The 
resulting E. coli lysate (approximately 260 mL) was placed into a 
dialysis bag (10 kD molecular weight cut-oil) in 4 L buffer (10 mM 
Tris-HCl, pH 7.4, containing 0.5 mM EDTA). The dialysis buffer 
was changed twice over a 16-24 h period. At the end of the 
dialysis, the cell debris was removed as a pellet by centrifugation at 
10000X g for 30 min at 10 °C. 

Dialysed and centrifuged recombinant /J 2 -AR supernatant was 
applied to a Ni + Sepharose 6 Fast Flow column that had been 
equilibrated with dialysis buffer (10 mM Tris-HCl, pH 7.4, 
containing 0.5 mM EDTA). The column was eluted with an 
increasing salt gradient of up to 500 mM imidazole in buffer A 
(20 mM potassium phosphate, pH 7.4, containing 500 mM NaCl) 
at a speed of 1 .0 mL/ min. The fraction of interest was collected by 
eluting the column using 50% buffer B (buffer A in the presence of 
500 mM imidazole, pH 7.4). Further purification of the collected 
fraction was performed on the Quaternary Sepharose Fast Flow 
column equilibrated with 18% buffer C (buffer A containing 
800 mM NaCl, pH 7.4). The /? 2 -AR solution was collected by 
gradient elution from the column using buffer C ranging from 
18%— 40% with a flow rate of 5.0 niL/min. The solution was 
finally desalted by the Sephadex G25 column. 

Oriented immobilisation of /? 2 -AR 

The /? 2 -AR was oriented immobilised on macroporous silica gel 
using a diazo coupling method (Figure 2). In detail, aminopropyl- 
bonded silica was synthesised according to the method in our 
previous report [28]. Two grams of aminopropyl-bonded silica gel 
was suspended in a /j-nitrobenzaldehyde methanol solution 
(18.2 mg/L, 100 mL) through stirring for 48 h at 60 °C. The 
resulting gel was filtrated by a sintered filter funnel and washed 
three times with methanol. The remnant methanol was removed 
by drying the gel in a vacuum oven at 60 °C for 12 h. The double 
bond in the produced gel was reduced to a single bond through 
reaction between the gel and 100 mL 8.6 |ig/mL NaBH 4 at room 
temperature for 12 h [29], using water as a solvent. After being 
filtrated and washed three times by water, the product was dried 
until the weight remained constant at 1 10 °C under vacuum. The 
nitro group in the gel was reduced to an amino group by mixing 
the gel with a SnCl 2 ethanol solution (8.1 mg/L, 100 mL) and 
reacting for 6 h at 90 °C. The suspension was cooled to 0 °C and 
alkalified with 50 mL 5% NaHC0 3 , followed by a sequential 
cleaning with saturated NaCl and distilled water [30]. The 
produced phenylamine group in the gel was then transformed to 
its diazonium form through reaction with excess HC1 (20 mM, 
30 mL), in which a freshly prepared NaN0 2 solution (20 mM) was 
added dropwise until reaching the turning point of potassium 
iodide starch paper changing to blue [31]. The whole procedure 
for preparing the diazonium salt was performed at 4 °C, and a 
post stirring for 10 min was needed to accomplish the reaction. 
The purified His-tagged /J 2 -AR was immobilised on the resulting 
gel by immediately adding the receptor (61.67 ng/mL, 10 mL) 
solution into the suspension and reacting for 2 h at 4 °C. Finally, 
the gel was washed three times with phosphate buffer (20 mM, 
pH 7.4). 

Chromatography 

The column containing the oriented immobilised /i 2 -AR was 
prepared by packing the stationary phase into a stainless steel 
column (2.1 mmx5.0 cm) using the slurry method at a pressure of 
4.0 xlO 7 Pa. Phosphate buffer (20 mM, pH 7.4) was used as the 
slurry solution. The control column containing exposed silica gel 



and attached histidine was prepared using the same synthesis 
method as that of the j8 2 -AR column. 

The mobile phases were prepared by adding 0 to 40 uM of 
ephedrine hydrochloride as the competing agent to 20 mM 
phosphate buffer (pH 7.4). All of the mobile phases were filtrated 
with 0.45 |im Millipore HV filters and degassed by ultrasonication 
before HPAC analysis. Working solutions of ephedrine and 
pseudoephedrine were prepared by dissolving the drugs in 
phosphate buffer (20 mM, pH 7.4) to a concentration of 
1.0 mM. The detection wavelength was 220 nm with an injection 
of 5.0 |jL for each drug. The flow rate was 0.2 mL/min for all of 
the analyses. Under the proposed chromatographic conditions, the 
void time was 1 .6 min using sodium nitrite as the injection solute. 
The capacity factor (K) was calculated by detecting the retention 
time in triplicate for each injected sample using the equation 
K = (tR— t())/t 0 . When the injection solute was pseudoephedrine, 
Eq. (1) was used to plot the curve of l/(£'-X) and [A], whereas Eq. 
(2) was used to achieve the plot between 1 / \fk' and [A] when 
ephedrine was used as the analyte. 

Thermodynamic study 

The thermodynamic investigation of ephedrine and pseudo- 
ephedrine on oriented immobilised /? 2 -AR was performed at 10, 
20, 30, 37 and 45 °C. Eq. (3) and (4) were used to calculate the free 
energy, enthalpy and entropy during the binding of the two drugs 
to the receptor. The semi-empirical rules by Ross et al. [26] were 
utilised to reveal the binding mechanism process of the two drugs 
to /? 2 -AR. 

Molecular docking 

The sequence of Sus scrofa j8 2 -AR (Accession Number: 
AF000134) was obtained from Genbank. The query sequence of 
the target receptor was searched find out the related protein 
structure to be used as a template by the BLAST (Basic Local 
Alignment Search Tool) program [32] against Protein Data Bank 
database. The three dimensional homology model was calculated 
using crystal structural coordinate of template on the basis 
alignment of target and template sequence of /J 2 -AR according 
to the procedures illustrated by Discovery Studio 2.5 (DS 2.5, 
Accerlrys Software Inc., San Diego, CA). Structural evaluation of 
the receptor was performed using two programs named PRO- 
CHECK and Verify 3D. Once the three dimensional model was 
generated, energy minimisation was carried out for the receptor 
before dockings. Structure of ephedrine and pseudoephedrine 
were constructed by DS 2.5. CDOCKER program, implemented 
in the software platform of DS, was used for the docking study. 
The pose cluster radius was set to 0.5 with a value of top hits as 10. 

Results and Discussion 

Purification of his-tagged /f 2 -AR 

The total protein in the desalted /? 2 -AR solution was determined 
by bicinchoninic acid (BCA). Using BSA as a standard, the 
calibration curve was plotted to assay the total protein in the 
solution. The linear regression equation was y = 0.0002x+0.2702, 
with a correlation coefficient of 0.9957. The concentration of the 
total protein in the /? 2 -AR solution was (6 1 .67 ±4.48) ug/mL. 
Sodium dodecyl sulphate polyacrylamide gel electrophoresis (SDS- 
PAGE) was used to detect the molecular weight of the protein in 
the desalted /? 2 -AR solution (Figure 3). The lanes of the desalted 
/? 2 -AR solution presented only one clear band corresponding to 
the band of the marker with a molecular weight of 66.5 kDa, 
based on which the molecular weight of the protein in the desalted 
/? 2 -AR solution is believed to be 66.5 kDa [30]. The purity of /? 2 - 
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Figure 2. Diagram for oriented attaching His-tagged /2-AR on the surface of diazonium salt coated silica gel. 

doi:1 0.1 371 /journal.pone.0094955.g002 



AR in the solution was calculated to be 95% of the total protein by 
densitometry quantification of the lane using the BandScan 
software. Furthermore, the desalted /J 2 -AR solution was also 
analysed by exclusion-size chromatography, which did not present 
any interference of the peak with a molecular weight of 66.5 kDa 
and gave a purity of 96% for the peak corresponding to /? 2 -AR. 
Taken together, the concentration of /? 2 -AR was approximately 
59 ug/mL. 

Immobilisation of His-tagged /? 2 -AR 

Random and oriented immobilised methods are the main assays 
for attaching a protein to solid substrates. Comparing with 
random methods, the proposed approach is oriented, and thus, it 
had the properties of better maintaining the bioactive character- 
istics and increasing the binding sites of the protein [33]. The 
concentration of the immobilised His-tagged /? 2 -AR on silica gel 
was calculated by detecting the free concentration of /? 2 -AR in the 




Figure 3. SDS-PAGE analysis of /2-AR solution. Lane 1, collection 
from Ni 2+ Sepharose 6 Fast Flow column eluted by 50% buffer B; Lane 2, 
protein molecular weight marker (the molecular weights from top to 
bottom are 1 16.0, 66.2, 45.0, 35.0, 25.0, 18.4 and 14.4 kDa, respectively); 
Lane 3-4, collection from Quaternary Sepharose Fast Flow column 
gradient eluted by 18%-40% buffer C; Line 5-9, collection from 
Quaternary Sepharose Fast Flow column eluted by 50%, 60%, 70%, 80% 
and 90% buffer C; Lane 10, /! 2 -AR supernatant. 
doi:1 0.1 371 /journal.pone.0094955.g003 



reaction solution before and after the diazonium reaction through 
the bicinchoninic acid method. The amount of yS 2 -AR attached on 
the gel was 200. 1 7± 7.47 pmol/g. This value was much greater 
than for the random immobilisation reported earlier [24,33]. The 
commonly reported oriented method for protein immobilisation 
included the use of metal chelate resin [34] and horseradish 
peroxidase [35]. This work developed a new oriented method 
through the reaction between the histidine tag in the protein and 
the diazonium salt activated silica gel. Compared with the 
previous two types of oriented methods, the proposed assay had 
more potency in varied protein immobilisation, attributed to the 
lack of structure selectivity of the protein. 

Specificity and stability of fe-AR column 

The specificity of the oriented /? 2 -AR column was tested by 
determining the retention times of specific drugs on the receptor. 
The retention times of the specific ligands, including salbutamol, 
clorprenaline, tobuterol and bambuterol were 5.0, 9.9, 12.5 and 
49.1 min, respectively, whereas the specific ligands of the /?r 
adrenoceptor and a [-adrenoceptor, esmolol and terazosin, both 
gave retention times of 1.6 min on the column (Figure 4A). The 
varied retention times illustrated that His-tagged /? 2 -AR still had 
the specificity for recognising its ligands after the immobilisation 
and packing procedures. Meanwhile, the rank order of the 
retention times of salbutamol, tulobuterol and clenbuterol is in line 
with the order of their association constants calculated by the 
radio-ligand binding assay [36]. These data confirms that it is 
reasonable to characterise the /? 2 -AR column using the retention 
times of the ligand targets of the receptor. 

The stability of jS 2 -AR was investigated by detecting the 
retention times of salbutamol, clorprenaline and tobuterol for 30 
consecutive days. The relative standard deviations (RSD) of the 
retention times on those days were 2.43%, 1.48% and 3.13%, 
indicating that the retention times of the drugs remained stable 
over the 30 days. The peak profiles were also stable over the 30 
days. In conclusion, the ;8 2 -AR column had good stability for at 
least 30 days. The specific drugs of /? 2 -AR give retention times 
approximate to the void time (Table 1) on the control column, 
which indicates that the retention of the specific ligands on the /? 2 - 
AR column is a consequence of the |8 2 -AR presentation but not a 
nonspecific interaction with the silica gel. 
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Figure 4. Representative chromatograms of control drugs, 
ephedrine and pseudoephedrine. A, Control drugs were eluted 
from /? 2 -AR column using 20 mM phosphate buffer (pH 7.4) as mobile 
phase at 0.2 mL/min. Non-retained peaks: sodium nitrite, terazosin, 
esmolol. Retained peaks: salbutamol, clorprenaline, tobuterol, bambu- 
terol (left to right); B, Ephedrine eluted from fe-AR column when the 
mobile phase containing ephedrine at different concentration (from 
right to left: 4, 8, 12, 20, 24, 30 and 34 (iM); C, Pseudoephedrine eluted 
from /)VAR column when the mobile phase containing ephedrine at 
different concentration (from right to left: 2, 4, 8, 12, 16, 24 and 30 \iM). 
doi:1 0.1 371 /journal.pone.0094955.g004 



Self competitive studies 

Self-competitive studies were performed where ephedrine 
served as both analyte and competing agent at 37 °C 
(Figure 4B). In this case, the injection volume of the analyte was 
5.0 uL of 40 uM ephedrine, and the competitive agent in the 
mobile phase was 4, 8, 12, 20, 24, 30 or 34 uM ephedrine. 
Triplicate injections were performed for each concentrations of the 
competitive agent in the mobile phase. The capacity factors of 
ephedrine in the /? 2 -AR column are summarised in Table 2. The 
capacity factors of ephedrine gave a negative response to the 
concentrations of ephedrine in the mobile phase. Meanwhile, the 
capacity factors did not reach saturated adsorption during the 
tested concentrations, providing a situation of the response in the 
linear scope of the isotherm curve. According to eq. (2), plotting 
the curve of 1 / \fk' against the concentrations of ephedrine in the 
mobile phase at 37 °C (Figure 5A), a good linear relationship is 
obtained with a regression equation of y = (2.58±0.02) X 10' x+ 
(0.43±0.01) and a correlation coefficient of 0.9962 over the seven 
data points, providing a proof of one type of binding site for 
ephedrine binding to ;8 2 -AR. Based on which ephedrine is believed 
to have only one type of binding site on the /J 2 -AR column. By 
dividing the intercept by the slope of Figure 5A, the association 
constant of ephedrine on the column was calculated to be 
(5.94±0.05)xl0 3 /M (Table 3). Through the slope of Figure 5A 
and the values of the association constant, the concentration of the 
binding site for ephedrine targeting /? 2 -AR was estimated to be 
(8.92±0.06) xlO" 4 M. 

Direct competitive study 

A direct competitive study was conducted using ephedrine as a 
mobile phase additive and pseudoephedrine as an injection solute 
at 37 °C (Figure 4C). As shown in Table 2, a decreasing capacity 
factor of pseudoephedrine was found when the concentration of 
ephedrine in the mobile phase was increased. The capacity factors 
of pseudoephedrine were not saturated during the tested 
concentrations, which indicated that the pseudoephedrine absorp- 
tion properties were in line with the Langmuir model. Eq. (1) was 
used to plot the curve between the pseudoephedrine capacity 
factors and the concentrations of ephedrine in the mobile phase. 
Using Eq. (1), a linear relationship was obtained between the 
l/(k'-X) of pseudoephedrine and the concentrations of ephedrine, 
as shown in Figure 5B. The value of X was calculated to be 
2.49±0.26 by non-linear regression. By applying the slope and 
intercept of the linear curve in Figure 5B, the association constant 
of ephedrine was determined to be (5.98±0.24)xl0 3 /M, as shown 
in Table 3. This value was in good agreement with the association 
constant obtained by the self- competitive studies. According to the 
hypothesis of zonal elution assay, the concentration of the binding 
sites for pseudoephedrine was (8.92±0.06) xlO + M based on the 
value of ephedrine calculated by self-competitive studies. Using the 
binding sites and the intercept of Figure 5B, the association 
constant for pseudoephedrine binding to /? 2 -AR was determined to 
be (3.80±0.02) xl0 3 /M, as shown in Table 3. 

The association constant of ephedrine is greater than that of 
pseudoephedrine even though the structures are similar. The 
difference between the structures of the two drugs is the 
configuration of the amino group, which produces the different 
spatial selectivity of the two drugs to the binding sites on /? 2 -AR 
during the interactions with the receptor. Although the association 
constants of the two drugs have the same order of magnitude, the 
exact values are varied due the high repeatability and precision of 
the HPAC assay, on which the spatial position of the amino groups 
are believed to contribute greatly to the binding of the two drugs 



PLOS ONE | www.plosone.org 



5 



April 2014 | Volume 9 | Issue 4 | e94955 



Binding of Ephedrine and Pseudoephedrine to ft-AR 



Table 1. 


Retention times of /? 2 -AR specific ligands 


on different columns. 






Llgand 


Retention times (min) 








Exposed silica gel 


Histidine coated gel 


Immobilized /' 2 -AR 


Salbutamol 


1.65 


1.62 


5.0 


Tobuterol 


1.58 


1.67 


9.9 


Clorprenaline 1.71 


1.64 


12.5 


Bambuterol 


1.56 


1.59 


49.1 


Notes: All the capacity factors are determined triplicate at 37 °C. 
doi:1 0.1 371 /joumal.pone.0094955.t001 


The void time is measured as 


1.6 min using sodium nitrite as an injection. 



and /? 2 -AR. This explanation can be confirmed by the research of 
Vansal et al. [37] who found that i^^-ephedrine has a maximal 
response relative to isoproterenol (EC50) of 72 |J.M, whereas 
7>?,2.S'-pseudoephedrine gives an EC50 value of 309 |tM on /? 2 -AR. 
Taken together, it is concluded that the varied affinity of the two 
drugs to /3 2 -AR may be one of the key factors for their different 
pharmacological potency through the signal pathway involving /? 2 - 
AR. 

Thermodynamic studies 

Thermodynamic studies between the two drugs and /? 2 -AR 
were performed at 10, 20, 30, 37 and 45 °C. Zonal elution was 
used to calculate the association constants of ephedrine and 
pseudoephedrine at those temperatures. According to eq. (3), 
linear relationships were found by plotting the IiiKa of ephedrine 
and pseudoephedrine against 1000/T, with correlation coefficients 
of 0.9716 and 0.9771 (Figure 6), respectively. The AH 8 and AS® 
for ephedrine binding to /J 2 -AR were calculated to be 
-(6.51±0.69) kj/mol and 50.94±0.31 J/mol-K. By eq. (4), the 
free energy during the binding was — (22.33±0.04) kj/mol 
(Table 4), indicating a spontaneous reaction during the binding 
procedure. For the interaction between ephedrine and /? 2 -AR, 
AH G <0, AS e >0 and AG 9 <0, which confirmed an exothermal 
process synergistically controlled by increasing entropy and 
decreasing enthalpy. During the process, AH 9 contributed 
29.15% to AG 9 , whereas AS 9 accounted for 70.85% of the 
contribution to AG 9 at 37 °C. Electrostatic interaction was 
believed to be the main force of the binding of ephedrine to /? 2 - 
AR. 

The values of AH e , AS 9 and AG 9 for the binding of 
pseudoephedrine to /? 2 -AR were —(7. 48 ±0.56) kj/mol, 
44.13±0.01 J/mol-K and — (2 1 . 1 7 ±0.02) kj/mol (Table 4). 
Similar to the interaction between ephedrine and /? 2 -AR, the 
binding of pseudoephedrine to the receptor was found to be a 
spontaneous reaction due to the negative value of AG 9 . The 
binding interaction proved to be mainly driven by electrostatic 
forces because AH 9 <0, AS 9 >0 and AG 9 <0 during the binding. 
The free energy for the binding of ephedrine and /? 2 -AR is more 
negative than that of pseudoephedrine, indicating that ephedrine 
more easily targets /? 2 -AR. This may be another reason for the fact 
that the two drugs have varied pharmacological potency. 

The association constants of ephedrine and pseudoephedrine 
binding to /? 2 -AR decreased with increasing temperature, whereas 
the binding sites presented a positive response to an increasing 
temperature. When the temperature changed from 10 °C to 45 
°C, the association constant of the two drugs decreased from 
7.17xl0 3 /M to 5.18xl0 3 /M and from 4.93xl0 3 /M to 
3.36xl0 3 /M, respectively. Meanwhile, the binding sites increased 
from 8.07 xl0~ 4 M to 9.48 xl0~ 4 M during the interaction, and 



the capacity factors of the two drugs on the /? 2 -AR column 
diminished from 5.40 to 4.63 and from 6.36 to 5.60 when the 
concentrations were 4.0 u\M. The association constants decreased 
1.4 and 1.5 fold, whereas the binding sites increased 1.2 fold. At 
the same time, the capacity factor diminished 1.2 fold. The 
retention behaviour of a solute on an affinity column ascribed to 
the synergistic contribution of the association constant and the 
binding sites. For ephedrine and pseudoephedrine binding to /? 2 - 
AR, the decreasing retention is mainly attributed to the 
diminishing association constants. The reason for this phenome- 
non may be related to the conformation change induced by the 
temperature. 

In the literature, it is reported that /? 2 -AR has three types of 
binding sites. One is the salt bridge between the amine group of 
the ligand and the carboxylate side chain of Asp 113 in the third 
hydrophobic domain. The other two types of binding interactions 
proceed through hydrogen bonds between the hydroxyl side chain 
of Ser 204 and the meta-hydroxyl group of the ligand or between the 
hydroxyl side chain of Ser 207 and the />ara-hydroxyl group of the 
ligand [38] . Ammonium salts mainly interact with the first type of 
binding site [39]. Being members of the group of ammonium salts, 
electrostatic interaction is believed to be the driving force during 
the binding of the two drugs to /? 2 -AR. In this presentation, both 
ephedrine and pseudoephedrine gave a thermodynamic profile of 
AG 9 <0, AH 8 <0 and AS 9 >0. This result confirmed that 
electrostatic interaction is the main force during the bindings. 
Thermodynamic studies in this work illustrated the same point to 
aforementioned report, which paves the way to the application of 
the thermodynamic method in revealing the binding mechanism 
of the interaction between ligands and the receptor. 

As noted in the literature [37], stereoselective and rank 
differences exist among the direct effects of ephedrine isomers 
on /J 2 -AR; li?,25-ephedrine is the most potent. In this case, the 
affinity of l-ff,2.S'-ephedrine binding to /? 2 -AR should be stronger 
than the other isomers. In this work, the association constants of 
l.ff,2.S'-ephedrine and lS,2£-pseudoephedrine were calculated as 
(5.94±0.05)xl0 3 M" 1 and (3.80±0.02) xlO 3 M"\ respectively. 
This result is notable ascribed to the high agreement with 
aforementioned report. Waldeck et al [39] has discovered that 
ephedrine has a marked positive inotropic effect on the papillary 
muscle with a pD 2 of about 5.95, which corresponds an 
dissociation constant of 1 1 200 nM for the drug binding to /? 2 - 
AR. Subsequent studies by Samama et al [40] and Seifert et al 
[41] have reported an dissociation constant of 8000 nM for 1^,26*- 
ephedrine binding to wild-type /? 2 -AR and a value of 
5100± 1300 nM to /J 2 -ARGsas. Based on the association constant 
in this work, the dissociation constant for lii',2.S'-ephedrme to the 
receptor was 168000 nM. Since the receptors used in these work 
are cloned with various gene sequence, it is expected that the exact 
values of dissociation constant for the drug differ from each other. 
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0.28 I ' 1 1 1 1 1 1 1 

0 8 16 24 32 

Ephedrine [uM] 

Figure 5. Change in 1 \ of ephedrine and 1 /(k'-X) of pseudoephedrine with mobile-phase concentrations of ephedrine. Data were 
from an oriented immobilised fe-AR column at 37 °C. The equations of the best fit lines were y = (2.58± 0.02)x10 3 x+(0.43 ±0.01) for ephedrine with a 
coefficient of 0.9962 and y = (1.76± 0.07)x10 3 x+(0.29 ±0.01) for pseudoephedrine with a coefficient of 0.9915. The value of X used for the plot of 
pseudoephedrine was set at (2.49±0.26). A, Ephedrine at different concentrations was added in the mobile phase. The injection solute was 5.0 nmol 
of ephedrine; B, Ephedrine at different concentrations was added in the mobile phase while 5.0 nmol of pseudoephedrine was injected as analytes. 
doi:1 0.1 371 /journal.pone.0094955.g005 



Molecular docking 

One template had been identified (PDB ID: 3KJ6, methylated 
/? 2 -AR-Fab complex) through homology searches in PDB by 
comparing and aligning the target sequence with that in PDB. The 
sequence identity was 50%. Followed by structural-based align- 
ment between query and target receptor, C-alpha atoms were 
selected for the model building ascribed to the highest identity 
(Figure 7 A). Under the desired modelling condition, ten annotated 
model structure of the receptor were predicted. The best model 
structure of /? 2 -AR superimpose to the template was believed to be 
pigAR.BLOOO 10001 (Figure 7B). It showed the lowest value in 



PDF (Probability Density Functions) Total Energy (-2866.93), 
PDF Physical Energy (-5004.22) and DOPE (Discrete Optimized 
Protein Energy) score (—31734.13). These results indicated that 
the best model of /i 2 -AR compared to the other models. The 
structural evaluation of /J 2 -AR model revealed that sterochemical 
and geometrical parameters incorporated in PROCHECK were 
satisfied in this model. The Ramachandran Plot of the receptor 
model calculated by PROCHECK was showed in Figure 7C. It 
was notable that most of the amino acid residues of /? 2 -AR 
distributed in the favoured region. Further verification of the 
model by Verify 3D confirmed that a good consistent of /? 2 -AR 
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Table 2. The 


capacity factors of ephedrine an 


d pseudoephedrine at 37 °C. 






Ephedrine (uM) 




Ephedrine (uM) 


k'p 


4 


5.01 ±0.02 


2 


5.87±0.01 


8 


4.84 ±0.02 


4 


5.81 ±0.01 


12 


4.58±0.03 


8 


5.73±0.01 


20 


4.31 ±0.01 


12 


5.66±0.01 


24 


4.08 ±0.01 


16 


5.54±0.03 


30 


3.80±0.01 


24 


5.47 ±0.03 


34 


3.64±0.01 


30 


5.38±0.01 



Notes: k'e and k'p represent the capacity factors of ephedrine and pseudoephedrine respectively, k'e is determined when ephedrine serves both injecting solute and 
competitive agents with the concentrations of 4, 8, 12, 20, 24, 30 and 34 ^M. While k'p is obtained by direct competitive study for which pseudoephedrine is used as an 
injection and ephedrine acts as competitive agents with the concentrations of 2, 4, 8, 12, 16, 24 and 30 uM. 
doi:1 0.1 371 /journal.pone.0094955.t002 



Table 3. The association constants of ephedrine and pseudoephedrine binding to /? 2 -AR at 37°C. 



Injection analyte 



Ephedrine 



Pseudoephedrine 



Ka ephedrine = (5.98±0.24) X 1 0 3 M 1 

Ka pseudoephedrine = 0.80 ±0.02) XlO 3 NT 1 



doi:1 0.1 371 /journal.pone.0094955.t003 



with the respective amino acid due to 65.30% of the residues 
having an averaged compatibility of an atomic model (3D) with its 
own amino acid sequence (ID). 

Site directed docking results showed that only one type of 
binding domain is engaged in the interaction between the two 



drugs and /? 2 -AR (Figure 7D). Asp 7 (corresponding to the 
reported Asp 1 13 of the reported /? 2 -AR) proved to be the exact site 
for the two isomers binding to the receptor. Meanwhile, no 
hydrogen bond was found during this interaction. It is widely 
known that Asp is an acidic amino acid with an isoelectric point of 



s 




1000/T 



Figure 6. Van't Hoff plots for the interaction of ephedrine and pseudoephedrine binding to^-AR. Circle, best-fit line with an equation 
of y = (0.78±0.08)x+(6.12±0.28) for ephedrine and a coefficient of 0.9716. Triangle, Best-fit curve for pseudoephedrine with an equation of 
y = (0.90±0.07)x+(5.31 ±0.23) and a coefficient of 0.9771. All the experiments were performed at 10, 20, 30, 37 and 45 °C. 
doi:1 0.1 371 /journal. pone.0094955.g006 
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Table 4. The free energy and association constants of ephedrine and pseudoephedrine in different temperature. 





Temperature (K) 


Ephedrine 




Pseudoephedrine 






Association constant 

(xio 3 wr 1 ) 


AG 8 (kJ/mol) 


Association constant 

(xio 3 wr 1 ) 


AG 9 (kJ/mol) 


283.15 


7.17±0.14 




4.94±0.08 




293.15 


6.69±0.10 


-(22.33 ±0.04) 


4.19±0.05 


-(21.17±0.02) 


303.15 


6.21 ±0.10 




4.02 ±0.05 




310.15 


5.94±0.05 




3.80 ±0.02 




318.15 


5.18±0.13 




3.35±0.07 





doi:1 0.1 371 /journal.pone.0094955.t004 



2.98. In this case, Asp possesses negative charge in a solution with 
pH of 7.4. While in the same environment, both ephedrine and 
pseudoephedrine exist as positive ions. Based on these points, we 
concluded that the binding of the two drugs to /? 2 -AR is driven by 



electrostatic interaction. This result is in good agreement with the 
results from thermodynamic investigation in this work and the 
report by Strader et al [42-43] . 




imiMiKimiM 



m 






mm 

mm 




Figure 7. Three dimensional model of Sus scrota /i-AR by homology modelling and its binding to ephedrine and pseudoephedrine. 

A, Structural-based alignment of /i 2 -AR (PDB ID: 3KJ6). The sequence shaded in cyan represents sequence similarity. Red colour means the 
transmembrane domain; B, The best predicted model of ^ 2 -AR (pigAR.BL00010001); C, Ramachandran plot of the /i 2 -AR model calculated by 
PROCHECK; D, Docking complexes overview of ephedrine//J 2 -AR and pseudoephedrine//j 2 -AR. Green illustrates Ephedrine. Blue illustrates 
Pseudoephedrine. 

doi:1 0.1 371 /journal.pone.0094955.g007 
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Conclusion 

A diazo coupling method was developed for oriented immobi- 
lising /? 2 -AR onto the surface of macroporous silica gel. Compared 
with random immobilised assays, the current method proves to 
have the advantages of better maintaining the bioactive charac- 
teristics of /? 2 -AR due to the greater amount of the immobilised 
protein on the gel and the increased binding sites on the surface of 
the receptor for recognising its ligands. The stationary phase 
prepared by the oriented method was also applied in exploring the 
interaction between jS 2 -AR and two drugs, ephedrine and 
pseudoephedrine. The varied association constants of the two 
drugs are believed to be the key reason for their different 
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